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Abstract: We study the nature of the finite-temperature chiral transition in QCD 
with Nf hght quarks in the adjoint representation (aQCD). Renormahzation-group 
arguments show that the transition can be continuous if a stable fixed point exists 
in the renormalization-group fiow of the corresponding three-dimensional theory 
with a complex 2Nf x 2Nf symmetric matrix field and symmetry-breaking pattern 
S\J{2Nf) — » S0(2A^j). This issue is investigated by exploiting two three-dimensional 
perturbative approaches, the massless minimal-subtraction scheme without e expan- 
sion and a massive scheme in which correlation functions are renormalized at zero 
momentum. We compute the renormalization-group functions in the two schemes to 
five and six loops respectively, and determine their large-order behavior. 
The analyses of the series show the presence of a stable three-dimensional fixed point 
characterized by the symmetry-breaking pattern SU(4) — >■ SO (4). This fixed point 
does not appear in an e-expansion analysis and therefore does not exist close to 
four dimensions. The finite-temperature chiral transition in two-fiavor aQCD can 
therefore be continuous; in this case its critical behavior is determined by this new 
SU(4)/SO(4) universality class. One-flavor aQCD may show a more complex phase 
diagram with two phase transitions. One of them, if continuous, should belong to 
the 0(3) vector universality class. 
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1. Introduction 

The thermodynamics of Quantum Chromodynamics (QCD) is characterized by a 
transition from a low-temperature hadronic phase, in which chiral symmetry is bro- 
ken, to a high-temperature phase with deconfined quarks and gluons (quark-gluon 
plasma), in which chiral symmetry is restored. See, e.g., Refs fl], |, ^, ^ for recent 
reviews. Although deconfinement and chiral symmetry restoration are apparently re- 
lated to different nonperturbative mechanisms, they seem to be somehow coupled in 
QCD. Indeed, lattice computations show that the Polyakov loop has a sharp increase 
around the critical temperature where the chiral condensate vanishes [Q. However, 
the interplay between the two effects is not clear yet. 

Insight into this question may be gained by investigating QCD-related mod- 
els, such as SU(Ac) gauge theories with Nj Dirac fermions in the adjoint repre- 
sentation (aQCD). In aQCD with Nf massless flavors the chiral symmetry group 
is S\J{2Nf), which is expected to be spontaneusly broken to S0(2A^j) at low tem- 
perature [|, 1^, H, 3, due to the presence of a nonzero quark condensate. In 
addition, unlike QCD with fermions in the fundamental representation, aQCD is 
also invariant under global Z^r^ transformations corresponding to the center of the 
SU(A'c) gauge group, as is the case in pure SU(Ac) gauge theories. This symmetry 
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breaks down in the high-temperature deconfined phase. Therefore, in aQCD one ex- 
pects two finite-temperature transitions: a deconfinement transition at Td associated 
with the breaking of the Ztv^, symmetry, and a chiral transition at character- 
ized by the symmetry-breaking pattern SU(2A^j) — >• S0(2iVj). Of course, it is also 
possible to have only one transition if the corresponding critical temperatures coin- 
cide. We should mention that this issue is interesting only for the values of Nf for 
which aQCD is asymptotically free. Since the first coefficient 60 of the (3 function 
[3{g) = —bog"^ + 0{g^) of SU(A''c) gauge theories with Nf adjoint Dirac fermions is 
60 = (11 -4A^/)A/'c/(247r2), aQCD is asymptotically free only for Nf < 11/4, i.e., for 
Nf = l and Nf = 2. 



Monte Carlo simulations of SU(2) [1^, |TT|] and SU(3) |]T2| gauge theories with 
two adjoint fermions show that the deconfinement and chiral transitions are well sep- 
arated with Trf < Tc. In the case of three-color aQCD, the Monte Carlo simulations 
reported in Ref. [1^ show actually that the ratio T^/Ta is rather large, Tc/T^ ~ 8, sug- 
gesting a rather weak interplay between the corresponding underlying mechanisms. 
They provide a rather clear evidence that the deconfinement transition associated 
with the center symmetry Z3 is of first order. Moreover, the available data at the chi- 
ral transition are apparently consistent with a continuous transition. However, they 
should be considered as rather preliminary and not conclusive, since a careful analy- 
sis of finite-size effects and of the approach to the continuum limit has not been done 
yet. The phase diagram of two-color aQCD in the temperature-chemical-potential 
plane has been recently discussed in Ref. [ p^ . 

In this paper we investigate the nature of the finite-temperature chiral transi- 
tion in aQCD, i.e. in four- dimensional SU(A^c) gauge theories with adjoint fermions, 
using renormalization-group (RG) arguments. Our study parallels the ones reported 
in Refs. |jl^, |T^, in which the nature of the finite-temperature transition in QCD 
with fermions in the fundamental representation was investigated. We consider ef- 
fective three-dimensional Landau-Ginzburg- Wilson (LGW) theories for the low- 
momentum critical modes associated with the bilinear quark condensate, which are 
described by a symmetric complex 2Nf x 2Nf matrix field, and look for stable fixed 
points (FPs) that may be associated with the symmetry-breaking pattern relevant 
for aQCD: SU{2Nf) S0{2Nf) [if the axial U(l) symmetry is effectively restored at 
Tc the symmetry-breaking pattern would be \J{2Nf) 0(2A^/)]. If such a stable FP 
does not exist the transition is of first order; otherwise, the transition may be con- 
tinuous or of first order, depending whether the system is or is not in the attraction 
domain of the stable FP. We study the RG fiow of the LGW theories in two field- 
theoretical (FT) perturbative approaches: the minimal-subtraction scheme without 
e expansion (in the following we will indicate it as 3(i-MS scheme) |0 and a mas- 



sive zero-momentum (MZM) renormalization scheme fl^. In the 3(i-MS scheme one 
considers the massless (critical) theory in dimensional regularization [|18[ , determines 
the RG functions from the divergences appearing in the perturbative expansion of 
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the correlation functions, and finally sets e = A — d = 1 without expanding in powers 
of e (this scheme therefore differs from the standard e expansion |19|). In the MZM 
scheme one considers instead the three-dimensional massive theory, corresponding to 
the disordered (high-temperature) phase, and determines the renormalization con- 
stants from zero-momentum correlation functions. We compute the /3 functions per- 
turbatively to five loops in the 3(i-MS scheme and to six loops in the MZM scheme. 
They are resummed by using a conformal-mapping method [^, 2T| taking into ac- 
count their large-order behavior, determined by means of the standard semiclassical 
analysis of instanton solutions. Comparison of the results of these two perturbative 
schemes provides a nontrivial check of the reliability of our conclusions. 

We briefly summarize our main results. In both 3(i-MS and MZM schemes the 
three-dimensional SU(4) LGW theory relevant for aQCD with two Dirac flavors 
shows a stable FP, that corresponds to a new three-dimensional universality class 
characterized by the symmetry-breaking pattern SU(4) S0(4). The correspond- 
ing critical exponents are rj ^ 0.2 and u ~ 1.1. Note that this FP does not appear in 
an e-expansion analysis (e = 4 — d) and therefore does not exist close to four dimen- 
sions. It is found only in genuinely three-dimensional analyses. This FT result implies 
that the finite-temperature chiral transition in two-flavor aQCD may be continuous. 
In this case, it belongs to the above-mentioned three-dimensional SU(4)/SO(4) uni- 
versality class. However, this does not exclude a first-order transition for systems 
that are outside the attraction domain of the stable FP. Note that, although SU(4) 
is locally isomorphic to SO (6), this universality class is definitely different from the 
vector 0(6) one, whose symmetry-breaking pattern is S0(6)— i>S0(5). No stable FP 
is found for the U(4) LGW $^ theory which should be relevant in the case the axial 
U(l) symmetry is effectively restored at T^, in this case, the transition would be of 
first order. For Nf = 1 the phase diagram may be more complex because the phase 
diagram of the corresponding SU(2)-symmetric effective theory has several transition 
lines joining at a multicritical point. In the parameter region relevant for aQCD, the 
possible phase diagrams have one or two phase transitions. One of them would be 
associated with the symmetry-breaking pattern SU(2)/Z2 = S0(3) S0(2). There- 
fore, if continuous, the transition would belong to the standard 0(3) (Heisenberg) 
universality class. The second transition, which could be absent in aQCD, would cor- 
respond to the symmetry-breaking pattern Z2 ® SO (2) — > SO (2), and, if continuous, 
would correspond to the Ising universality class. Note that the symmetry breaking 
SU(2) — > S0(2) occurs in this case through two different phase transitions. 

The predicted critical behavior for two-flavor aQCD should be compared with 
that of two-flavor QCD with quarks in the fundamental representation. Also in 
this last case the finite-temperature chiral transition may be continuous, although 
it would belong to a different universality class, the vector 0(4) universality class. 
Monte Carlo simulations of lattice QCD |^2|, ^ ^, ^ seem to be consistent with 
a continuous transition in the 0(4) universality class, although they are not yet 
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sufficiently precise to be conclusive. 

The paper is organized as follows. In Sec. |^ we derive the LGW $^ theory 
relevant for the finite-temperature chiral transition in aQCD, using universality and 
RG arguments. In Sec. |^ we investigate the RG flow of the LGW theory with U{N) 
symmetry, which would be relevant for aQCD if the U(l) anomaly were effectively 
suppressed at Tc. We report our perturbative calculations in the 3(i-MS and MZM 
renormalization schemes and their analyses. In Sec. ^ we consider a more general 
LGW theory in which the U(A^) symmetry is explicitly broken to SU(A^). We discuss 
the one-flavor case and present a perturbative analysis for = 4, which is relevant 
for two-flavor aQCD. The appendix contains an analysis of the vacuum structure of 
the LGW (j)'^ theories relevant for aQCD. 

2. The effective Landau- Ginzburg- Wilson model at the chiral 
transition 

In the vanishing quark-mass limit, the fermionic part of the QCD Lagrangian with 
Nf adjoint Dirac fermions is given hj C = tpf'-ff^D^ipf, where D^^^ = 9^5'^''+zy4^/"* 
and are the generators of the adjoint representation, i.e. the structure constants. 
Using the antisymmetry of the structure costants, one may rewrite the Lagrangian 
in terms of two-component Weyl spinors as C = i lyJcx^D^Wj, where = 
(Id, a). Id is the identity matrix, and cTj are the Pauli matrices; see, e.g., Refs. 
[0, |, §. The actual symmetry is U{2Nf) = \]{l)A^[S\J{2Nf)/Z2Nf], which is larger 
than the symmetry U(A^/)ij U(A^/)l of QCD with fermions in the fundamental 
representation. The U(1)a subgroup is anomalous at the quantum level and thus the 
symmetry reduces to S'U{2Nf). 

At T = the symmetry is expected to be spontaneously broken due to a nonzero 
quark condensate {ipip). As a consequence of the Pauli principle, the quark bilin- 
ear condensate must belong to the symmetric second-rank tensor representation of 
S\J{2Nf), which has dimension 2Nj + Nf. Condensation along one of its directions 
gives rise to the symmetry breaking 

SU(2Ar^) S0(2Ar^) (2.1) 

and to 2A^| + Nf — 1 Goldstone modes. See, e.g., Refs. H, ^ for more details. 
With increasing the temperature, a phase transition characterized by the restoring 
of the chiral symmetry is expected at a given T^; above Tc the quark condensate 
vanishes. Therefore, the finite-temperature phase transition is characterized by the 
symmetry-breaking pattern ( |2.1| ) and a complex symmetric 2A^/ x 2Nf matrix order 
parameter In the case the U(l) symmetry is effectively restored at Tc, the 

relevant symmetry-breaking pattern would be 

U(2A^^) ^ 0(2Ar^). (2.2) 
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This possibility is however rather unhkely. Indeed, instanton calculations in the 
high-temperature phase ||2^ suggest that the axial U(l) symmetry is not restored at 
Tc, in analogy with what happens in lattice QCD with fermions in the fundamental 
representation ||2^, and as also suggested by the finite-temperature behavior of the 
topological susceptibility in the pure SU(Ai'c) gauge theories, see, e.g., Ref. p6| . 

In order to investigate the nature of the finite-temperature transition in aQCD 
with Nf light fiavors, we follow the reasoning already applied in Refs. [H, ^ to 
the study of the finite-temperature transition in QCD with light fermions in the 
fundamental representation. 



(i) Let us first assume that the phase transition at Tc is continuous (second or- 
der) for vanishing quark masses. In this case the critical behavior should be 
described by an effective three-dimensional {3-d) theory. Indeed, the length 
scale of the critical modes diverges approaching Tc, becoming eventually much 
larger than 1/Tc, which is the size of the euclidean "temporal" dimension at 
Tc. Therefore, the asymptotic critical behavior must be associated with a 3-d 
universality class characterized by a complex symmetric 2Nf x 2Nf matrix or- 
der parameter $jj and by the symmetry-breaking pattern (|2.1|) [or ( p. 21) if the 
U(l) symmetry is effectively restored at Tc]. 

(ii) According to RG theory, the existence of such universality classes can be in- 
vestigated by considering the most general LGW $^ theory for a complex 
symmetric N x N matrix field $jj with = 2Nf and the desired symme- 
try and symmetry-breaking pattern. The most general U(A^)-symmetric LGW 
Lagrangian containing up to quartic terms in the potential is 

= Tr (9^<l>t)(9^$) + rTr $t$ + ^ (Tr + ^Tr (<l>t$)' . (2.3) 

Stability requires uq + vq > and Nuq + vq > 0. The symmetry group of 
this Lagrangian is U(A^). First, it is invariant under the transformations $ 
U^U^ , where [/ is a unitary U(A^) matrix, i.e., under the group U{N)/Z2 
(the quotient Z2 is due to the fact that matrices ±f/ give rise to the same 
transformation). Second, it is invariant under $ — $^ The symmetry group 
is therefore Z2®U(A^)/Z2 = U(A^). For vq > 0, Lagrangian ( p.3| ) shows the 
expected vacuum structure and symmetry-breaking pattern, see appendix. In 
the low-temperature phase the potential is minimized by taking 

^-,.oc(«'^^). (2.4) 

where Id is the A^/2-dimensional identity matrix. Note that $min is invari- 
ant under vector U{Nf)v transformations, in agreement with the Vafa-Witten 
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theorem ||^. The symmetry of the vacuum is O(A^) and thus the symmetry 
breaking pattern is U(A^) 0{N). 

The axial anomaly reduces the symmetry to SU(A^) and thus new terms must 
be added. The most relevant one is proportional to (det $^ + det $) , which is 
a polynomial of order N in the field $. For > 4 such a term is irrelevant at 
the transition. Instead, for N = 2 and = 4, the determinant must be added 
to Lagrangian ( |2.3| ), obtaining 

jOsu{n) = CuiN) + Wo (det $t + det $) . (2.5) 

For N = 2 additional terms depending on the determinant of $ should be 
added, in order to include all terms with at most four fields compatible with 
the symmetry. They will be discussed in Sec. [4.1|. The symmetry of La- 



grangian ( |2.5| ) is Z2 ® SU(A^). Indeed, it is invariant under transformations 
$ U^U'^ , with U unitary and det U = ±1. Taking into account that U and 
—U correspond to the same transformation, the invariance group is SU(A^). 
Moreover, the model is invariant under the Z2 transformations $ ^ $^ For 
^0 > — 3|wo|/2, the vacuum structure is identical to that discussed for the U(A^) 
theory, satisfying the Vafa-Witten theorem 0. The corresponding symmetry- 
breaking pattern is Z2 ® SU(A^) — >■ Z2 ® SO(A^). Note that the symmetry 
group contains an additional Z2 with respect to Eq. ( |2.1| ). This additional 
invariance, related to the transformation $ — is a consequence of the her- 
miticity (reality) of the effective Lagrangian and corresponds to the discrete 
parity symmetry of the aQCD Lagrangian. This additional Z2 is not broken 
at the transition, so that the relevant symmetry-breaking pattern is indeed 
SU(A^) SO(A^), as discussed before. It must be noted that, if a 6* term is 
added to the aQCD Lagrangian, this additional symmetry is not present and 
the effective Lagrangian is no longer hermitian. 

(iii) The critical behavior at a continuous transition is described by the stable FP 
of the theory, which determines the universality class. The absence of a stable 
FP indicates that the phase transition is not continuous. Therefore, a nec- 
essary condition of consistency with the initial hypothesis that the transition 
is continuous, cf. (i), is the existence of stable FPs in the theories described 
by Lagrangians (p.3|) and ( p.5|) . If a FP exists, the transition may be either 
continuous, belonging to the universality class associated with the stable FP, 
or of first-order, if the system is not in the attraction domain of the FP. 

The effective theories ( p.3|) and (|2.5|) are defined for any A^. However, for aQCD 
they are relevant only for even and A^ < 4, i.e., for N = 2 and A^ = 4. Indeed, 
aQCD is asymptotically free only for Nj < 11/4, i.e., for A^ < 11/2. 
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3. Renormalization-group flow of the U(A^) LGW theory 



3.1 Perturbative expansions 

In order to investigate the RG flow of Lagrangians ( |2.3| ) and ( p.5| ), we employ two 
different perturbative approaches: the MZM renormahzation scheme [|r^ and the 
3c/-MS scheme [|l6l; see Refs. [pi] , ^ for recent reviews. In the flrst case, one consid- 
ers the three-dimensional massive theory corresponding to the disordered phase, and 
expresses the zero-momentum renormahzation constants in terms of zero-momentum 
renormalized quartic couplings. In the second case, one considers the massless (criti- 
cal) theory in dimensional regularization within the minimal-subtraction scheme 
RG functions are obtained in terms of the renormalized couplings and of e = 4 — d. 
Subsequently e is set to its physical value e = 1, providing a three-dimensional scheme 
in which the 3-d RG functions are expanded in powers of the MS renormalized quar- 



tic couplings. This scheme differs from the standard e expansion |]T9| in which one 
expands the RG functions in powers of e. 

In order to renormalize the U(A^)-invariant theory ( p. 3D in the MS scheme, one 

sets 



(3.1) 



where u and v are the MS renormalized quartic couplings. The renormahzation con- 
stants Zfj,, Zu, and Z^ are normalized so that Z(i){u,v) ~ 1, Zu{u,v) ~ u, and 
Zy{u,v) ^ V at tree level. Here is a d-dependent constant given by A^ = 
2rf-i^rf/2p^^^2). Moreover, one deflnes a mass renormahzation constant Zt{u,v) by 
requiring ZtF^^'^-* to be flnite when expressed in terms of u and v. Here F^^'^-* is 
the one-particle irreducible two-point function with one insertion of Tr$^$. The /? 
functions are computed from 



du 



Uo,Vo 



dv 
(9/i 



Uo,Vo 



They have a simple dependence on d: 



Pu = {d- 4)u + Bu{u, v), (3y = {d- A)v + B^{u, v), 



(3.2) 



(3.3) 



where the functions Bu{u, v) and By{u, v) are independent of d. The FPs of the theory 
are given by the common zeroes of the (3 functions. Their stability is controlled by 
the eigenvalues of the matrix 



dpu/du dpu/dv 
d(3.u/du d(3.u/dv 



(3.4) 
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A FP is stable if all the eigenvalues of its stability matrix have positive real part. 
The RG functions t]^ and rjt associated with the critical exponents are defined by 

d In Z^^t 



(3.5) 



MO, DO 

The functions rjfp^t are independent of d. The standard critical exponents 77 and u are 
related to the RG functions rj^^t and to the location u*, v* of the FP by 

r] = r]^{u\ V*), iy=[2 + 7]t{u\ v*) - r]^{u\ v*)]-' ■ (3-6) 

We computed the MS series to five loops. For this purpose we used a symbolic 
manipulation program that generates the diagrams (305 up to five loops) and com- 
putes their symmetry and group factors, and the compilation of Feynman integrals 
of Ref. EH]. The functions Bu,v{u,v) defined in Eq. (|3.3|) are given by 



^ / X A^^ + A^ + S 2 , 3 2 9Ar2 + 9Ar + 42 , , 

Bu{u,v) = + {N + l)uv + -y 3.7 

4 4 lo 



ll(Ar + l) 2 5Ar2 + 15Ar + 92 , 3(Ar + 2) 3 „ , 

4 32 8 ^ 



u( ^ ^ 2^ + 5 2 5Ar2 + 5Ar + 82 2 ,^ 

BJu.v) = 6UV ^ V uv (3.8) 

4 16 

HAT + 20 2 3Ar2 + 21Ar + 60 3 sr^ , i , 

4 32 ^3 

The coefficients a^- and a^^- for 4 < i + j < 6, are reported in Tables |I] and || 
respectively. In order to save space, we report them numerically, although we have 
their exact expressions in terms of fractions and C, functions with integer argument. 
We do not report the series of the RG functions rj^^t, since they will not be used in 
the following. They are available on request. 

In the MZM scheme the theory is renormalized by introducing a set of zero- 
momentum conditions for the one-particle irreducible two-point and four-point cor- 
relation functions: 

rll,,,,,(p) = '^'L = '^'^i^.'^-''.^*"' + + o{p')] , (3.9) 

d9'aia20^bib2 



.(4) ... 



' (0) 

aia2,bib2,ciC2,did2 ^ ■' 



S'^lia2SH^b2^^crC2S'^did2 



(3.10) 



zero mom. 



2 4^ 327r 

= m g _|_ ^ _|_ ^2 ('^^oia2,bib2,ciC2,did2 ~^ '^^aia2,bib2,c\C2,dld2) i 

where F is the generator of the one-particle irreducible correlation functions (effective 
action), and 

Uaia2Mb2,ciC2,did2 2 i^aiCi^bidi^a2C2^b2d2 ~^ ^aidi^biCi^a2d2^b2C2) 5 
Vaj^a2,bib2,ciC2,did2 ~ 2 {^aiCi^bidi^a2d2^b2C2 ~^ ^aidi^biCi^a2C2^b2d2) ■ (3-11) 
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In addition, one introduces the renormalization function Zt{u,v) that is defined by 
the relation T^^^f^^ ,^^^^(0) = ^a^b^^aj^jZ^"^, where r*^^'^)(0) is the zero-momentum one- 
particle irreducible two-point function with one insertion of Tr$^$. The MZM f3 
functions are defined by 



f3u{u,v) = m 



du 
dm 



Pv{u, v) = m 



dv 
dm 



(3.12) 



Uo,Vo 



The RG functions r]^ and rjt associated with the critical exponents are defined by 

d In Z^^t 



d\nm 



(3.13) 



Mo, Do 



We computed the MZM RG functions to six loops. In this case we used the 
compilation of Feynman integrals of Ref. to compute the needed 1438 Feynman 
diagrams. The f3 functions are given by 



Puiu,v) 



2 4:{N+1 



-uv + 



-V + 



12 2N + b 



i+j>3 



uv + 



i+j>3 



(3.14) 
(3.15) 



The coefficients 6"^- and V^^ up to six loops, i.e. for 3 < i + j < 7, are reported in 
Tables ^ and ^ respectively. Here we do not report the six-loop series of the RG 
functions rj^^t] they are available on request. 

Our calculations can be checked by considering some particular cases. For = 1 
the two quartic terms are identical and one recovers the 0(2) vector model, while for 
f = one obtains the 0(A^^ + A^)-symmetric model. In these cases we can compare 
our perturbative expansions with those reported in Refs. ^ ^]: we find full 
agreement. In addition, for N = 2 the model is equivalent to an 0(2)®0(3) model 
Indeed, if one sets 



$ = ^(^fei + i^fc2)efc. 



(3.16) 



fc=i 



where (61,62,63) = i(Id, zo"i, zaa). Id is the 2x2 identity matrix, o"j are the Pauli 
matrices, and is a 3 x 2 real matrix, Lagrangian ( p.3| ) for N = 2 can be written 
as 



C 



91,0 
' 4! ' 



2 ^2 I 92fi_ 

4! 



2 \2 



(3.17) 
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^'^^ 3 3 3 

^1,0 = + -^vo, g2,o = --vq. (3.18) 

This model has an exphcit symmetry [0(2)®0(3)]/Z2 = U(l) ® SU(2)/Z2 = U(2), 
as expected. It has aheady been studied because it should describe transitions in 
frustrated spin models with noncollinear ordering, the superfiuid transition of ^He, 
etc., see, e.g., Refs. |2^ and references therein. Six-loop series in the MZM 
expansion |^ and five-loop series in the MS scheme |^ have already been computed 



for generic 0(2)®0(M) symmetric models. These perturbative series agree with 
ours, once we rewrite them in terms of the renormalized couplings corresponding to 



the bare ones defined in Eq. ( |3.18D . In the MS scheme, the relation is gi = |f , 



g2 = —ft", where gi and g2 are the MS renormalized coupling of model ( |3.17| ). In 



the MZM scheme, the correspondence is cji = u + v/2 and g2 = —dv/lA, where gi^2 
are the MZM couplings normalized as in Ref. |3^ (actually, there they were called u 
and v). 

3.2 Large-order behavior 

Since perturbative expansions are divergent, resummation methods must be used 
to obtain meaningful results. Given a generic quantity S{u, v) with perturbative 
expansion S{u,v) = "^ij CijU''v^ , we consider 



S{xu,xv) = '^^Sk{u,v)x'^, (3.19) 

k 



which must be evaluated at x = 1. The resummation of the series can be done by 
exploiting the knowledge of the large-order behavior of the coefficients, generically 
given by 

Skiu, v) ~ k\ [-A{u, v)f [1 + 0{k-^)] . (3.20) 

The quantity A{u,v) is related to the singularity tg of the Borel transform B{t) that 
is nearest to the origin: ts = ~1/A{u,v). The series is Borel summable for a; > 
if B{t) does not have singularities on the positive real axis, and, in particular, if 
A{u,v) > 0. A semiclassical analysis based on the instanton solutions, see, e.g., 
Ref. PT|, indicates that the MS and MZM expansions are Borel summable when 



In this region we have 



where 



u + v>0, Nu + v>0. (3.21) 



A{u, v) = c Max (u + v,u + (3.22) 



c=l (MS), 

1-32997 ^^-^^^ 
c= (MZM), 
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respectively for the MS and MZM expansions. Under the additional assumption that 
the Borel-transform singularities lie only in the negative axis, the conformal-mapping 
method outlined in Refs. I^O], |21| turns the original expansion into a convergent one 
in the region ( 3.21| ). Alternatively, one may use the Pade-Borel method, employing 



Fade approximants to analytically extend the Borel transform. 

Outside the region ( p. 211 ) the expansion is not Borel summable. However, for 
f > 0, if the condition u + v > holds, a naive extension of the results obtained for 
the Borel-summable case — but this a quite nonrigorous procedure — indicates that 
the Borel-transform singularity closest to the origin should still be on the negative 
axis. Therefore, the large-order behavior should still be given by Eq. ( p.20|) with 



A{u,v) given by Eq. ( p.22| ). Thus, by using A{u,v) as given by Eq. ( |3.22| ) and the 
conformal-mapping method, one may still take into account the leading large-order 
behavior, and therefore hope to get an asymptotic expansion with a milder behavior, 
which may still provide reliable results. 

We should mention that the MS series are essentially four- dimensional, so that 
they may be affected by renormalons that make the expansion non-Borel summable 
for any u and v, and are not detected by the semiclassical analysis leading to 
Eqs. ( |3.20|) , (|3.21| ) and (|3.22|) ; see, e.g., Ref. The same problem should also 



affect the MS series of 0{N) models. However, the good agreement between the 
results obtained from the FT analyses [0 and those obtained by other methods 



indicates that renormalon effects are either very small or absent (note that, as shown 
in Ref. |3^, this may occur in some renormalization schemes). For example, the 



analysis of the five-loop perturbative 3(i-MS series |T^ gives u = 0.629(5) for the 
Ising model and u = 0.667(5) for the XY model, that are in good agreement with the 
most precise estimates obtained by lattice techniques, such as u = 0.63012(16) |^ 
and z/ = 0.63020(12) || for the Ising model, and z/ = 0.67155(27) ^ for the XY 
universality class. On the basis of these results, we will assume renormalon effects 
to be negligible in our analyses of the 3c?-MS series. 

3.3 Analysis of the series 

One can easily identify two FFs in the theory described by the Lagrangian Cu(n), 
without performing any calculation. The first one is the Gaussian FF with u = v = 0, 
which is always unstable. Since, for vq = 0, Cu{n) is equivalent to the Lagrangian of 
the 0(A^^ + A^)-symmetric vector model, there is an 0(A^^ + A^) FF with v = and 



M > 0. The results of Ref. |^ on the stability of the 0(M)-symmetric FF under 
generic perturbations can be used to prove that, for any N >2, the 0(A^^ -|- A^) FF 
is unstable in theories ( p.3|) and ( p.5|) . Indeed, the term Tr ^ in the Lagrangian 

Cu{N), which acts as a perturbation at the 0(A^^ + A^) FF, is a particular combination 
of quartic operators transforming as the spin-0 and spin-4 representations of the 
0(A^^ + A^) group, and any spin-4 quartic perturbation is relevant at the three- 
dimensional 0(M) FF for M > 3 since its RG dimension 1/4 4 is positive for 
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M > 3. In particular, y^^i fa 0.27 at the 0(6) FP |42|, increases monotonically with 
increasing M and approaches the value 1/4,4 = 1 in the large-M limit. 

In order to investigate if other FPs are present, one may use the e expansion. 
A simple analysis, that requires only the one-loop terms in the expansion of the (3 
functions, indicates that for e ^ 1 only the Gaussian and 0(A^^ + N) FPs exist and 
that none of them is stable, in agreement with the preceding analysis. Thus, close 
to four dimensions, the transition is of first order. However, there may be FPs that 
exist in three dimensions but are absent for e <^ 1. This is indeed what happens 
in the Ginzburg-Landau model of superconductors, in which a complex scalar field 
couples to a gauge field and in the 0(2)®0(iV) theory describing the critical 
behavior of frustrated spin models with noncollinear order Ref. [p5| , and refer- 
ences therein, and the superfluid transition in ^He |^5[|. Thus, to correctly identify 
the three-dimensional critical behavior, we must employ strictly three-dimensional 
perturbative schemes. Therefore, we study below the RG flow by using the 3(i-MS 
scheme and the MZM scheme. 

For N = 2 Lagrangian ( p.3|) is equivalent to Lagrangian ( |3.17|) written in terms 
of a real 3x2 matrix field. This model has already been studied by FT methods 

H, |§] in the MZM (to six loops) and 3d-MS (to five loops) schemes. These 
studies found evidence of two stable FPs: one (called chiral FP) with attraction 
domain in the region (72,0 > [ P5| , M and another one (called collinear FP) in the 
region g2fi < [45|. According to the mapping ( 3.18| ), the domain f > relevant for 
aQCD corresponds to (72,0 < 0, and thus the collinear FP is the relevant one. Using 
the results of Ref. [0, we find therefore a stable FP at u* = —0.54(8), v* = 1.14(6) 
in the 3d-ME scheme and at u* = -3.0(3), v* = 4.6(3) in the MZM scheme.^ These 
results are confirmed by a direct analysis of the perturbative series in u and v. 

Let us now consider = 4 and investigate the possible existence of a stable 
FP with f > 0. Let us first consider the 3(i-MS scheme. In order to find the 
zeroes of the /5-functions, we first resum the expansions of Bu{u,v) and B^{u,v) 



^ We mention that a mapping similar to that reported in Eqs. (3.16) and ( |3.17D exists also for 
Lagrangian ( |2.3D when <I> is a generic (not necessarily symmetric) complex 2x2 matrix. Such a 
model has a larger symmetry group [U(2)i (g) U(2)fl;]/U(l)y and is relevant for two-flavor QCD 
fl^ , when the effect of U(l)yi anomaly is neglected. Setting <i> = J^^^i'^ni'^ni + ^^^'2), where 
CT^ = (Id, (Ti) and '5 is a 4x2 real matrix, Lagrangian (2.3) corresponds to Lagrangian ( 3.17 ) with 
9i,o — l^^o + jVq, 32,0 — "i'^o- The analysis of the perturbative series in the 3d-MS and MZM 
schemes to five and six loops respectively |Q shows that theory (3.17) for N — 4 has a stable FP 
with 52 < 0. Therefore, the 3-d [JJ{2)l(E>V{2)r]/V{1)v model has a stable FP, located at u* « -0.5, 
V* « 1.2 in the 3d-MS scheme and in u* = -3.4(3) and v* = 5.3(3) in the MZM scheme (we use 
here the normalizations of Ref. jl^). Note that this stable FP is not found close to four dimensions 
by analyses based on the e expansion ||lj, Q . Its existence implies that systems characterized by 
the symmetry-breaking pattern [U(2)i ® U(2)ij]/U(l)y — > Uy(2)/Uy(l) can undergo continuous 
transitions if they are in the attraction domain of this stable FP. This FP was overlooked in Ref. jl5| . 
Indeed, the analysis of the MZM expansions was limited to the region —2 ^ u,v < A, that does not 
include the FP reported above. 
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u 



Figure 1: Zeroes of the Sd-MS /5- functions for = 4 in the region v > (sohd hne 
for Pu, dashed hne for /3^,). The black dots indicate the Gaussian and the 0(20) FPs. 
The dotted hne 4u + u = separates the Borel-summable region (at the right) from the 
non-Borel-summable region (at the left). 



defined in Eq. ( p.3D . More precisely, we consider the functions Ru,y{u,v, x) = 
Bu,v{ux, vx) I x^ . For each function i?„ „ we consider several approximants constructed 
using the conformal-mapping and Pade-Borel methods. In particular, in the case of 
the conformal-mapping method, we consider different values of the resummation 



parameters a and h (typically a = —1, 0, 1, 2 and b = 2, . . . , 18), see Ref. for 
definitions. 

The zeroes of f3u and f3y in the region v > are shown in Fig. |I| for = 4. They 
were obtained by using the conformal-mapping method and parameters a = and 
6 = 6. Approximants corresponding to different values of a and b, or obtained by 
using the Pade-Borel method, give similar results. We find no evidence of additional 
FPs beside the Gaussian and the 0(20) ones in the v = axis. 

The analysis of the six-loop MZM series was done in a similar way. We considered 
the functions Pu{u, v)/u and /3„(m, v)/v and we apply the conformal-mapping method 
as before. The results are perfectly consistent with the Sci-MS ones. No additional 
FP is found for = 4. 

According to our FT results, a consistent model for a three-dimensional continu- 
ous transition characterized by a complex symmetric N x N matrix order parameter 
and symmetry-breaking pattern U(A^) 0{N) does not exist for = 4. Therefore, 
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the phase transition in such systems must be of first order. For N = 2 instead, 
there is a stable FP; thus, systems characterized by the symmetry breaking U(2)— 
0(2) can undergo a continuous transition, if they are in the attraction domain of the 
stable FP. 



4. Renormalization-group flow of the SU(A^) LGW theory 

We now study the effect of the breaking of \J{1)a- As we mentioned in Sec. ^ we 
must add to the U(A^) Lagrangian a term proportional to det $ + det $^ Such a 
term is irrelevant for > 4: in this case the critical behavior is described by the 
U(A^) theory. Therefore, we only consider the cases N = 2 and = 4, which are the 
only ones of interest for aQCD. Indeed, > 4 corresponds to Nf > 2 and aQCD is 
asymptotically free only for Nf < 11/4 = 2.75. 

4.1 The case N = 2 

In the case N = 2 {Nf = 1), the determinant is quadratic in the fundamental field 
$, and other terms must be added to the effective Lagrangian, i.e. all possible 
interactions with at most four fields. This leads to the LGW Lagrangian 

Csu{2) = Tr(a^$T)(a^$) +rTr$^<l> + Wo (det $^ + det $) + ^ (Tr<l>t$)^ 
+^Tr (<l>t$)' + ^ (Tr (det + det $) + j [(det ^^)^ + (det $)'] ,(4.1) 

where $jj is a complex symmetric 2x2 matrix. Note that there is no need to include 
also a term proportional to (det $)(det $^). Indeed, for any two dimensional matrix 
M, we have 2detM = (TrM)^ - TrM^. Therefore, 

2(det$)(det$1') = (Tr <l>^<l>)^ - Tr ($^$)^ (4.2) 

Since the Lagrangian ( [4.1| ) has two quadratic mass terms, it describes several tran- 



sition lines with a multicritical point. 

Model ( |4.1| ) can be written in terms of two real three-component fields (pi and ipi 
by writing $ = J2'j=ii^j + where (61,62,63) = (Id/2, 'jcTi/2, 2(73/2), Id is the 

2x2 identity matrix, and cTj are the Pauli matrices. One obtains 

+^{<pr + ^^{^r + |;0V + ^(0 ■ (4.3) 

with = r + WQ, = r-wo, ao = |(2uo + fo + 2a;o + i/o), bo = l{2uo + vo-2xo + yo), 



a model is studied in Ref . pDl . One finds a multicritical point with several transition 



Co = |(2uo + 3fo — Ho), and rfo = — 3(fo + yo)- The mean-field phase diagram of such 
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Figure 2: Possible phase diagrams for model ([4.3D . In phase diagrams (al), (a2), and (a3) 
the multicritical point corresponds to a second-order transition, while in phase diagrams 
(bl), (b2), and (b3) the multicritical transition is of first order. Thick lines represent first- 
order transition, while thin lines correspond to second-order transitions. We label with 
H, XY, and Is the lines corresponding to transitions belonging to the Heisenberg {N = 3 
vector), XY (N = 2 vector), and Ising universality classes, respectively. 

lines. The analysis reported in Ref. indicates that the multicritical point may 
either belong to the chiral 0(2) ® 0(3) universality class or may correspond to a 
first-order transition. Note that in Ref. |^ the full RG flow was not studied and 
thus one cannot exclude the presence of other stable FPs, beside the 0(2) (g) 0(3) 
chiral FP. This possibility is rather unlikely, but an explicit calculation of the RG 
flow in the full parameter space is needed to settle the question. Here, we will not 
consider this additional possibihty, since it is of no relevance for aQCD, as we discuss 
below. In the plane t, g, where g parametrizes the breaking of the symmetry (p (p, 
the possible phase diagrams are reported in Fig. |^. 

In the case of one-flavor aQCD, the U(2) theory obtained by setting wq = 0, 
Xq = 0, Ho = should correspond to the multicritical point because of the larger 
symmetry group. Therefore, the relevant FP at the multicritical point is the one 
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found for the U(2) theory in Sec. 3^. The other FPs of the theory (O) are of no 
relevance since they are not present in the U(2) model. Thus, if chiral symmetry is 
not restored at the transition, the behavior of aQCD as a function of temperature T 
corresponds to the behavior observed along a nontrivial line in the {t, g) plane in one 
of the phase diagrams reported in Fig. ^ ^ However, we do not know which of the 
phase diagrams applies. Generically, we expect two phase transitions. One of them, 
that should occur at higher values of T, may be either of first order or belong to the 
Heisenberg universality class [whose symmetry-breaking pattern is SO(3)^SO(2)], 
which has been accurately studied in the literature, see, e.g., Refs. [^, and 
references therein. Such a transition would be associated with the symmetry breaking 
Z2(S)SU(2) — > Z2®0(2). The lower-temperature transition, that may not necessarily 
exist, may be of first order, or continuous. According to Fig. 0, the transition may 
belong either to the Ising or to the XY universality class. It is easy to see that 
in the latter case the U(l)v symmetry would be broken, violating the Vafa-Witten 
theorem |Q. Therefore, the only possibility is an Ising transition that corresponds to 
the breaking Z2 Cg> 0(2) ^ Z2 C?) S0(2). Note that the symmetry breaking SU(2) 
S0(2) is realized here through two different transitions and that an SU(2)/SO(2) 
universality class does not exist. One may discriminate between the one- or two- 
transition hypothesis by determining the symmetry of the T = phase, which differs 
by a Z2 group in the two cases. 



4.2 The case = 4 

For = 4 {Nf = 2) the determinant is a quartic-order term, giving rise to a 
generalized LGW theory ( p. 5] ) with three quartic parameters. In this case stability 
requires 

Uo + Vo>0, Auo + Vo-2\wo\>0. (4.4) 

As discussed in the Appendix, the symmetry-breaking pattern and vacuum structure 
appropriate for aQCD are realized for vq > — ||wo|. Note that the FPs already 
identified in the U(4) theory, i.e. the Gaussian and the 0(20) FP, are also FPs of the 
SU(4) theory. They are both unstable and thus are of no relevance for the critical 
behavior. 

similar scenario with a phase diagram characterized by a multicritical point applies to the 
physically more interesting case of QCD with two flavors in the fundamental representation. On the 
basis of the results mentioned in footnote |l|, the presence of a stable FP in the [\J(2)l(^\J{2)ii]/U{1)v 
<i>^ theory (which was overlooked in Ref. has a direct consequence on the possible phase 

diagrams in the T,g plane, where g parametrizes the effective breaking of the \J{1)a symmetry. 
Indeed, it leaves open the possibility that at the multicritical point the transition is continuous, 
with an 0(4) critical line starting from it. In this case, even an infinitesimal breaking of the U(l)yi 
symmetry can give rise to an 0(4) critical behavior. This should partially correct the conclusive 
remarks of Ref. on the phase diagram relevant for two-flavor QCD. 
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In order to study the RG flow of tlie tlieory (|2.5|) for = 4, we consider tlie 
3(i-MS sclieme and the MZM scheme as before. In the 3(i-MS scheme we set 

$= [Z^{u,v,w)Y'^^R, (4.5) 
uo = Ad^i^Zu{u,v,w), 
Vq = AdfJ,^Z^{u,v,w), 

Wo = —AdiJ.'Z^{u,v,w), 

and determine the corresponding (3 functions to five loops. They are given by 

Pu = —eu + 7u + 5uv + -V + 3w u -u v -uv (4.6) 

4 8 4 4 ^^ 

i+j+k>i 

Py = —ev + 3uv + -j-v^ — 3w'^ — -^u'^v — IGuv'^ — 6v^ (4.7) 
4 8 



9 



+9uw + -VW + > 



cLi.u'v^w'', 



2 " / ^ 

«+i+fc>4 

3 91 2 7 2 7 

= —ew + 6UW VW u w V w uvw (4. 

^ 2 8 4 4 ^ 



i+j+A;>4 



The coefficients Cijk up to five loops, i.e. for 4<'i+j + A;<6, are reported in 
Table |. We also computed the five-loop expansion of the RG functions rj^^t defined 
as in Eq. (|3.5|): 

m= E (4-9) 

i+i+fc>l 

The coefficients ef-^. and e*^-;. up to five loops, i.e. for 1 < i + j + A; < 5, are reported 
in Table ^. The standard critical exponents are related to rj^^t by 

T] = T]4,{u\v\w*), u =\2 + iqt{u\v\w*)-r]^{u\v\w*)Y^ , (4.10) 

where m*, f *, w* are the coordinates of the stable FP in the case it exists. 

In the MZM scheme, beside the renormalization conditions (|3.9| ) and (|3.10|) , we 
also set 

.(4) 



aia2,bib2,ciC2,did2\ y 



^'^aia2 6<!>b^b2S<!>c^C2S<^d^ 



d2 



zero mom. 



2 4 d ^ 

~ -^Wea^biCidiea2b2C2d2, (4-11) 



-17- 



where eijki is the completely antisymmetric tensor (ei234 = 1). We computed six-loop 
series in the MZM scheme. The MZM f3 functions are given by 

P^ = -u + u^ + ^uv + + ^w^+ Y '^Iku'v'w^ (4.12) 

j+j+fc>3 
i+j+fc>3 

3 3 

(3^ = -w + -uw - —vw + Y dfj^u'v^w''. (4.14) 

i+j+k>3 

The coefficients dijk up to six loops, i.e. for 3 < i + j ' + /c < 7, are reported in 
Table 0. We also computed the RG functions ri(f,^t defined in Eq. ( p.l3| ) to six loops. 
The coefficients efjf^ and e*^;, defined in Eq. ( [4. 91) are reported in Table |. 

Again, close to four dimensions there are only two FPs, the Gaussian and the 
0(20) FPs (the latter located at -u = e/7, v = 0, and w = 0), i.e. those already 
found in the U(4) theory, which are both unstable. In order to investigate the 
possible existence of other FPs in three dimensions, we analyze the 3(i-MS and MZM 
perturbative expansions. Note that the model is invariant under the transformations 
$ — > 6*'^/''$, {u,v,w) — > {u,v,—w). This implies that Pu, Pv, and Pw/u) are even 
in w. Therefore, we can restrict our search of FPs to the w > space; if a FP 
with coordinates u,v,w > exists, there is also another FP with the same critical 
properties at u,v, —w. 



We use the conformal-mapping method already employed in Sec. p.3| and the 
large-order behavior of the perturbative series. Writing 

S{xu,xv,xw) = 'Y^Sn{u,v,w)x'^, (4-15) 

n 

semiclassical calculations based on instanton solutions give 

s„ ~n![-A(M,^;,w)]"n^ (4.16) 

V \w\ 



A{u, V, w) = cMax \ u + v, u -\ h 



where the constant c is given by c = 3/4 and c = 0.0474989 respectively for the MS 
and MZM expansions. The perturbative expansions are Borel summable for 

u + v>0, 4u + v-2\w\>0. (4.17) 

The 3(i-MS five- loop series are analyzed as in the U(A^) case and indicate the 
presence of a new FP with w ^ 0. Most approximants of the five- loop series (more 
than 90% of the approximants with a = —1,0,1,2 and b = 3, . . . , 18) present a 
common zero with w > 0. They are approximately 50% at four loops. This new 
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FP is located at u* = 0.23(9), v* = 0.34(6), and w* = 0.54(8), where the error 
takes into account the spread of the resuhs of the approximants considered. Note 
that the error (spread of the results) is rather large, essentially because this FP lies 
relatively far from the origin (note that the 0(20) FP lies at u* ~ 0.17), and also 
because it is close to the boundary of the Borel summable region. The analysis 
of the corresponding stability matrix shows that this FP is stable, and therefore it 
determines the universal properties of continuous transitions in systems described 
by the SU(4) LGW $^ theory (|2.5| ). We also computed the corresponding critical 
exponents by evaluating the RG functions 77^ 4 at the FP. From the analysis of the 
expansions of t]^ and (2 + r/^ — ?7<^)~^, we obtain r] = 0.23(8) and z/ = 1.1(3), where 
the error takes into account the spread of the approximants and the uncertainty 
on the FP coordinates. The presence of the stable FP is confirmed by an analysis 
based on Pade-Borel approximants; indeed, using [4/1] Pade approximants for all 
/3 functions, one obtains u* ^ 0.29(5), v* ^ 0.23(1), and w* ^ 0.51(4) (the errors 
indicate how the results change when the parameter b is varied between 4 and 18), 
which is substantially consistent with the results of the conformal-mapping analysis. 
Apparently, the Pade-Borel method gives a more accurate estimate of the location 
of the FP. But, note that the error is only related to the dependence on 6 of a single 
Pade approximant; thus, it is likely to be underestimated. As we have already seen 
in many other instances, the conformal-method error should provide a more realistic 
estimate of the real accuracy of the result. 

The rather low precision of the 3(i-MS results may give rise to some doubts on 
the existence itself of the stable FP, and therefore it calls for an independent and 
nontrivial crosscheck. This is provided by the analysis of the perturbative expansions 
in the alternative three-dimensional MZM scheme. We follow the same steps as in 
the U(A^) case, considering the conformal-mapping method and approximants with 
a = —1,0, 1,2 and 6 = 3, ...18 (see Ref. for definitions). The analysis confirms 
the presence of a stable FP at u* = 0.0(3), v* = 3.5(2) and w* = 4.6(2), where the 
error is related to the spread of the results. The presence of a FP is stable with 
respect to the number of loops. Indeed, at five loops (resp. four loops) we find a 
FP at u* = 0.4(2), V* = 3.2(2) and w* = 4.3(2) (resp. u* = 0.7(3), v* = 3.6(5) and 
w* = 4.6(6)). For comparison, in this scheme the 0(20) FP is located at u* ~ 1.18 
[0. We also estimate the critical exponents, finding r] = 0.20(3) and z/ = 1.1(3) 
(again from the analysis of the expansions of r]^ and (2 + 77^ — 77^)"^), in substantial 
agreement with the results of the 3c?-MS scheme. Also in this scheme results are not 
very precise, although they are apparently more accurate than the 3(i-MS ones. 

In conclusion, the analysis of the perturbative expansions in the 3(i-MS and MZM 
schemes provides evidence for the existence of a stable FP, and therefore of a new 
universality class that describes continuous transitions in systems described by the 
LGW Lagrangian ( |2.5|) with = 4 and symmetry-breaking pattern SU(4) SO (4). 
Although the results of the analyses do not appear particularly precise, the substan- 
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tial agreement between the two schemes makes us confident on their rehabihty. Note 
that the existence of a FP does not exclude the possibihty of observing first-order 
transitions; indeed, this is still possible for systems outside the attraction domain of 
the stable FP. 

We recall that in Sec. p.3| we found no stable FP in the U(4) LGW $^ theory, 
which should be relevant if the U(l) symmetry broken by the anomaly is effectively 
restored at Tc, thus suggesting first-order transitions for the corresponding systems. 
The existence of a stable FP with w ^ in the SU(4) theory is an interesting example 
of the so-called phenomenon of softening of first-order transitions: by breaking some 
symmetry of the original model, a first-order transition may become a second-order 
one. This phenomenon is well known in spin systems: the introduction of quenched 
disorder may soften the first-order transitions of pure systems In disordered 



systems translational invariance is the broken symmetry. In our case, it is instead 
the internal symmetry of the system that gets reduced. 
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A. Vacua of the V{N) and SU(4) theories 

Let us consider the potential 

= rTr$t$ + ^ (Tr$t$)' + ^Tr(<l>^$)', (A.l) 

and let us determine the fields that minimize it. The stationarity condition gives the 
equation 

r^ab + y$a;,Tr <l>t$ + ^ (<f$^<f)^, = 0. (A.2) 
If $ is a solution of Eq. ( |A.2| ), it is easy to show that 

^f/(7V)($)L„i = ^Tr$t$. (A.3) 

In order to determine the solutions of Eq. (|A.2|) , let us perform a polar decomposition 
of $: we write $ = PU, where P = V <I>^$ is a hermitian positive semidefinite 
matrix and f/ is a unitary matrix. Then we diagonalize P, writing P = VPdV^ with 
Pd = diag(Ai, A2, . . . , Aiv), Ai > A2 > . . . > A^r > 0, and V unitary. Now, if k is the 
rank of $, we have for 1 < a < k the equation 

- + yTrPj + |A^ = 0, (A.4) 
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that shows that Xa does not depend on a. Summing this equation over 1 < a < k, 
we obtain 

Tr$t$ = TrP| = V = (A.5) 

a=l 

Since Uq + Vq > and Nuq + fo > 0, for r > we must have k = 0, i.e. the only- 
possible solution is $ = 0: for r > the system is disordered. For r < instead, any 
k with < k < N is acceptable. Using Eq. (|A.3| ) we obtain 



which shows that the energy depends only on the rank k of the solution. A simple 

calculation shows that the minimum of the potential is attained for A; = if f o > 
and = 1 if fo < 0. 
For k = N, we have 

2r 

P = -T-r Id, (A.7) 

where Id is the identity matrix. Thus, ^ oc Q, where Q is a symmetric unitary 
matrix. Of course, being U(A^) connected, we can write Q = [/2 = UU^ , with 
U unitary and symmetric. Thus, modulo symmetry transformations, we may take 

"^min ^ Id. 

U k = 1, = VaVb, whcrc Va is an A^-dimensional complex vector. Since, for 
any v there exists U G U(A^) such that Va = UabWb, where w = (a, 0, . . . , 0), with 
a real and positive, we can write $ as 

2r 

$ = UhU^, (A.8) 

Mo + Vq 

with Ji = diag(l, 0, . . . , 0). Thus, modulo symmetry transformations, we may take 

$mm oc Ii. 

Let us now determine the symmetry-breaking patterns. The theory is invariant 
under the transformations $ —>■ U^U^ , U unitary, which form a group isomorphic 
to U(A^)/Z2 (the quotient is due to the fact that U and —U give rise to the same 
transformation). The transformations that leave invariant $min oc Id are those with 
U G O(A^), so that the invariance group is 0(A^)/Z2. If $min oc Ii, the transforma- 
tions that leave invariant have the form U = (±l)©i?, R G U(A^— 1), so that the 
invariance group is [U(A^ - 1) Zsl/Zg = U(A^ - 1). Beside $ U^U'^, the model. 
Id, and Ii are invariant under the Z2 transformations Thus, the symmetry 

breaking pattern is V{N) O(A^) for vq > and V{N) U(A^- 1) OZa for vo < 0. 
Note that for A^ = 2 they are identical since 0(2) = S0(2) O Z2 = U(l) ® Zg. 

To conclude, let us note that, as a consequence of the Vafa-Witten theorem 
[^], the vector group U{Nf)y can never be broken. Therefore, the vacuum must be 
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invariant under the U{N/2) transformations 

with U eU {N/2). It is easy to verify that^ the only symmetric matrices Q such that 
UvQUy = Q are proportional to 

where Id is the N/2x N/2 identity matrix. Thus, the only possible vacuum for aQCD 
is proportional to J. It is easy to verify that this is possible for f > (as we showed 
any matrix proportional to a symmetric unitary matrix is a possible vacuum), but 
not for V < 0. Thus, in the aQCD case we must restrict ourselves to the case v > 0. 
Now let us consider the effect of the determinant for = 4. The potential is 

Vsu(N){^) = Vu(N){^) + wo{det^ + det^^). (A.12) 

The stationarity condition gives then 

Cldiei 02^2^2 

(A.13) 

Now, we write again $ = VPdViU = VPd[{detVy/^V^UV*]V^ with Pd diagonal 
and V — y(dety)~^/^. Thus, modulo a symmetry transformation, we can simply 
write $ = PdU, with U unitary. Substituting in the previous equation we obtain the 
equation 

rXa + ^XaTr Pj + ^A^ - «;o(det C/) J] = 0. (A.14) 

c^a 

The term proportional to is the product of three eigenvalues and thus it is relevant 
only if the rank of $ is 3 or 4. If the rank is three, assume A4 = 0. The previous 
equation for a = 4 gives 

A1A2A3 = 0, (A.15) 
that contradicts the hypothesis A1A2A3 7^ 0. Thus, there is no solution with rank 

^Indeed, if we write 

(A.10) 

where a, b, and c are N/2 x N/2 matrices, we must have Uall^ = a, UcU^ = c, UhU'^ = b. The 
first two conditions imply a = c = (it is enough to consider infinitesimal matrices U), while the 
third one gives b oc Id. 



(A.9) 



(A.ll) 
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Assume now that the rank k is 4. Since all Xa are real and nonvanishing, (det U) 
must be real, hence equal to s = ±1, since U is unitary. Let us now determine the 
eigenvalues. Eq. (|A.14| ) implies 



Rab = XaEa — Xb^b — (A^ — A^) 



r + 



Mo 



^0, 



2TrP| + y(A^ + A^ 



0. 



(A.16) 



Now, consider R12 = 0, -R13 = 0, and -R23 = 0. It is immediate to verify that these 
equations imply that at least two eigenvalues among Ai, A2, A3 are equal. Without 
loss of generality we assume Ai = A2. An analogous discussion indicates that at 
least two eigenvalues among Ai, A3, A4 are equal. Thus, there are two possible cases: 
Ai = A2 = A3; Ai = A2 and A3 = A4. Finally, consider i?24 = 0; we obtain that either 
A2 = A4 or 2r + MoTr P| + vo{Xl + A|) = 0, that gives a relation between A2 and A4. 
In conclusion, we obtain two different classes of solutions: 



(i) all eigenvalues are equal with 
and potential 



-2r 



4mo + + 2swo 



sol 



4mo + vo + 2swo 
(ii) one eigenvalue differs from the others: 



(A.17) 



(A.18) 



A? — Xl 



2swf] 



A 



-Ai, 



-2rMo 



3mom5 + v^ + 4moWo + 4f 



2 ' 



(A.19) 



with potential 



+ Aw^ 



sol 



SMo-Wq + ""o + 4MoWo + '^VoW^ 



2 • 



(A.20) 



Note that the solution with s = 1 exists only if WqVq is positive, the one with 
s = — 1 in the opposite case. For wq = this solution corresponds to the rank-3 
solution we have determined before. 



Comparing the value of the potential for the different solutions, we obtain finally: 

(a) for mq < and |m7o| < — 2fo/3 the relevant solution has rank k = 1 and thus, 
modulo symmetry transformations, we can take $min oc Ji. The symmetry of 
the vacuum is SU(3)®Z2. Such a case is not of interest for aQCD, since it is 
not invariant under U(2)y. 
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(b) for wq > and wo > —2vo/3 the relevant solution is given by Eq. ( A. 17 ) with 



potential ( |A.18| ), setting s = — 1. Therefore, U is symmetric with det U = —1. 
It follows that U = e^'^I^W'^ ^ V G SU(4). Modulo symmetry transformations, 
one can therefore choose $min e^'^^Hd. The vacuum is invariant under SO (4) 
(note that the symmetry $ ^ is broken here). This case is not relevant for 
aQCD since the vacuum breaks $ ^ 



for Wq < and |wo| > —'^vq/S the relevant solution is given in Eq. (IATtI) with 



potential ( |A.18|) , setting s = 1. Therefore, U is symmetric with detU = 1. It 
follows that U = VV'^ with V G SU(4). Modulo symmetry transformations, 
one can therefore choose $min oc Id. The vacuum is invariant under 80(4)® Z2. 
The matrix J is one of the possible vacuum solutions and thus this case is of 
relevance for aQCD. 
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4,0 


10.7397 + 2.927707V + 


2.992167V2 + 


0.1289067V3 + 


0.06445317V4 


3,1 


17.5092+ 18.12647V + 


1.234377V2 + 


0.6171877V3 




2,2 


25.2599 + 9.506987V + 


3.952527V2 + 


0.01171877V3- 


f 0.002929687V4 


1,3 


13.4014 + 8.45267Ar + 


0.6269537V2 - 


f 0.1171877V3 




0,4 


3.61949+ 1.546277V + 


0.3011697V2 






5,0 


-63.5743 - 21.06067V 


- 22.36387V2 


- 2.605927V3 - 


- 1.301947V4 




+0.001220707V5 + 0.000406901Ar^ 






4,1 


-140.774- 153.3877V 


- 25.23527V2 


- 12.64257V3 - 


- 0.02983217V4 




-0.009944047V5 








3,2 


-260.744- 143.3887V 


- 66.57727V2 


- 2.011647V3 - 


- 0.526260A/"' 


2,3 


-212.910- 157.2627V 


- 23.22537V2 


- 5.144077V3- 


- 0.02237617V4 




-0.003442487V5 








1,4 


-107.874 - 60.89917V 


- 15.91067V2 


- 0.7090647V3 


- 0.07845487V4 


0,5 


-20.6853 - 12.88277V 


- 2.21799Ar2 


- 0.2159297V3 




6,0 


475.552 + 183.0117V + 200.15267V2 - 


f 34.55617V3 + 


17.96137V4 + 0.8182707V^ 



+0.2699067V6 - 0.002442877V^ - 0.000610718Ar^ 

5.1 1339.84 + 1524.757V + 373.2697V2 + 195.2517V3 + 10.32007V4 + 3.40600iV5 
-0.02915457V6 - 0.007288647V^ 

4.2 3025.04 + 2080.71A^ + 1065.77iV2 + 97.11367V3 + 29.09287V4 + 0.178866iV5 
+0.02834667V6 

3.3 3355.06 + 2824.447V + 654.5687V2 + 172.5647V3 + 4.471047V4 + 0.769792Af5 

2.4 2471.77 + 1715.817V + 548.5347V2 + 48.23677V3 + 6.775867V4 + 0.02001497V5 
+0.001956717V6 

1.5 944.177 + 692.2877V + 168.5637V2 + 25.20917V3 + 0.7603727V4 + 0.0526262Ar5 
0,6 157.642 + 109.2787V + 29.09447V2 + 2.788407V3 + 0.1637817V4 



Table 1: 3d-MS scheme: coefficients cf. Eq. {^). 
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4.0 

3.1 25.4497 + 2.33904iV + 2.288iV2 - 0.1015627V3 _ 0.0507812Ar4 

2.2 39.1784 + 23.3634iV - 0.0214843iV2 - 0.136718iV3 

1.3 29.5574+ 12.9326iV + 2.32341A^2 

0,4 7.12024 + 3.75368A^ + 0.486716A^2 ^ o.0253906Af3 



5,0 







4,1 


-177.097 - 22.30327V - 22.5027iV2 _ o.392090iV' 


^- 0.1785297V4 




+0.0210195A^5 + G.00700652iV6 




3,2 


-374.999 - 241.790iV - 14.5218iV2 - 4.728817V3 


+ 0.1790727V4 + 0.0584779A^5 


2,3 


-422.833 - 232.645iV - 55.3180iV2 + 0.1763229iV3 o.l375007V4 


1,4 


-209.997 - 132.747iV - 23.9377A^2 _ 2.37050iV3 




0,5 


-42.8890 - 25.88087V - 5.72264iV2 - 0.396930iV' 


^ - 0.01360947V4 


6,0 







5,1 


1513.47 + 243.234iV + 250.177iV2 + 13.8095iV3 + 


6.714367V4 




-0.228540A^5 _ o.07627337V6 - 0.00007995647V'^ 


- 0.00001998917V8 


4,2 


4108.25 + 2830.727V + 327.3697V2 + 121.9347V3 - 


0.9550587V4 




-0.471074A^5 _ 0.005697307V6 - 0.001090387V7 




3,3 


6155.99 + 3929.91A^ + 1110.53A^2 ^ ss.assSA^^ _,. 


4.607237V4 




+0.0693464A^5 _^ o.00888516iV6 




2,4 


4661.51 + 3367.087V + 766.2607V2 + 98.83457V3 - 


0.03282127V4 - 0.02475457V5 


1,5 


1909.14 + 1333.31A^ + 360.926A^2 _^ 37.io827V3 + 


2.380667V4 


0,6 


317.910 + 228.597A^ + 59.13417V2 + 7.38414Af3 ^ 


0.3185977V4 + 0.007276537V5 



Table 2: 3d-MS scheme: coefficients cf. Eq. (|!|). 
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{N^ + N + 




3,0 


-28.1481 - 


6.07407 iV- 


6.07407 




2,1 


-29.6296 - 


29.6296 iV 






1,2 


-30.8148- 


S.llllliV- 


1.7037 




0,3 


-8-4A^ 









4.0 199.64 + 54.9404 N + 56.2893 + 2.69789 + 1.34894 

3.1 329.228 + 342.521 N + 26.5875 iV^ + 13.2938 

2.2 476.927 + 191.216 N + 81.4672 A^^ _,_ 1.24378 + 0.310945 

1.3 254.359 + 165.976 N + 16.1265 + 3.27131 

0,4 68.2248 + 29.746 + 5.95053 A^^ 

'Jfi -1832.21 - 602.521 N - 638.341 - 71.4848 - 35.3533 iV^ + 0.466938 
+0.155646 

4.1 -4128.34- 4492.49 iV- 725.59 Ar2_ 

3.2 -7702.2 -4371.43 TV- 1991.42 iV2- 

2.3 -6312.85 -4744.76 TV- 712.149 iV2 

1.4 -3175.43- 1760.63 iV- 458.141 7V2 
0,5 -588.999- 356.825 iV- 57.9832 



356.03 + 8.11729 + 2.70576 
31.3703 iV3- 3.61396 iV* 

- 143.849 - 2.7972 iV* - 0.430339 N" 

- 26.2356 7V3_ 2.65955 iV* 

- 4.8188 TV^ 



6,0 

5,1 

4,2 

3,3 
2,4 

1,5 
0,6 



20770.2 + 7819.56 N + 8488.12 + 1340.35 A^^ 
+3.54529 iV6 + 0.204945 iV^ + 0.0512362 iV^ 
59871.8 + 67528.7 TV + 15357.4 iV^ + 7788.99 
+5.00417iV6 + 1.25104 iV^ 
136493. + 94965.3 N + 46606.9 iV^ + 2966.99 
+9.60758 iV^ 
151865. + 128949. TV - 
111193. + 76504.5 TV - 
+0.068185 iV6 
41631.6 + 30068.3 TV - 
6838.45 + 4597.45 TV- 



29147.1 TV2 
24548.1 TV2 

7403.69 TV2 
1180.76 TV2 



7117.85 TV3 
2500.48 TV3 

1147.35 TV^ 
105.622 TV3 



678.319 TV4 + 9.91857 TV^ 

135.84 TV'i + 51.1183 TV^ 

844.792 TV4 + 37.7957 TV^ 

261.112 TV4 + 54.5075 TV^ 
355.428 TV4 + 2.33913 TV^ 

43.3788 TV4 + 2.1912 TV'^ 
5.48048 TV^ 



7,0 
6,1 
5,2 
4,3 
3,4 
2,5 
1,6 
0,7 



-271300. - 114181. TV - 126851. TV^ - 25604.3 TV^ - 13465.7 TV^ - 793.197 TV^ 
-259.288 + 4.52143 TV'^ + 1.30604 TV** + 0.117121 TV^ + 0.0234242 TV^" 
-956030.- 1.1157 10^ TV- 323182. TV2- 171156. TV^- 11419.1 TV* - 3707.12 TV^ 
+89.4175 TV^ + 27.8058 TV^ + 3.63425 TV« + 0.72685 TV^ 

-2.57728 10^ - 2.06957 10^ TV - 1.09407 10^ TV^ - 129858. TV^ - 38311.9 + 107.385 TV^ 
+97.7575 TV^ + 29.1378 TV^ + 6.06393 TV« 

-3.6309 10^ - 3.40291 10^ TV - 1.00539 10*^ TV^ - 271968. TV^ - 9603.86 TV* - 1333.52 TV^ 
+ 103.325 TV6 + 23.2754 TV^ 

-3.48896 10^ - 2.80102 10^ TV - 1.03092 10^ TV^ - 132759. TV^ - 19077.3 TV* - 122.197 TV^ 
+ 17.666 TV^ - 5.06322 10"^ TV^ - 4.70902 10"* TV** 

-1.96674 10^ - 1.61821 10^ TV - 481214. TV^ - 88493. TV^ - 5662.61 TV* - 448.949 TV^ 
-8.28982 TV^ - 0.500033 TV^ 

-641804. - 491982. TV - 154530. TV^ - 21660.9 TV^ - 2089.08 TV* - 70.0898 TV^ 
-3.10875 TV^ 

-89165.7 - 68401.1 TV - 19266.8 TV^ - 2704.06 TV'^ - 149.873 TV* - 5.22722 TV^ 



Table 3: MZM scheme: coefficients 5"-, cf. Eq. (|3l^ ) 
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{N^ + N + 




3,0 











2,1 


-54.8148- 


3.40741 iV- 


3.40741 




1,2 


-53.6296- 


29.6296 iV 






0,3 


-20.1481- 


7.111117V- 


1.03704 





4.0 

3.1 469.334 + 41.8539 N + 40.6028 - 2.50221 - 1.25111 TV* 

2.2 720.915 + 426.443 N - 4.51405 - 4.50686 

1.3 544.203 + 235.122 TV + 40.6104 iV2 

0,4 131.416 + 69.5514 N + 9.46842 iV^ + 0.560179 
~dfi 

4.1 -5032.69 - 584.288 iV - 584.021 - 0.0404379 - 1.45685 TV* - 1.72396 
-0.574653 iV^ 

3.2 -10619.9 - 6724.27 iV - 257.525 TV^ _ 80.9869 - 7.06796 iV* - 3.10306 iV^ 

2.3 -12009.2 - 6419.67 iV - 1446.44 TV^ + 3.57388 TV^ _ 4.26467 TV* 

1.4 -5985.6 - 3722.89 iV - 685.886 - 78.9515 iV^ 

0,5 -1243.21 - 746.976 iV- 168.232 7V2 - 11.7457 TV^ _ 0.496927 A'"' 
~6fi 

5.1 64749.3 + 9324.6 TV + 9463.87 7V2 + 274.9 TV^^, 128.057 iV*- 12.1619 7V^ 
-5.53845 iV^^- 1.27242 iV'^- 0.318104 7V« 

4.2 175059. + 117375. N + 9362.48 iV^ + 3195.12 - 187.625 - 64.2153 
-8.29114iV6- 2.37504 iV^ 

3.3 263538. + 160906. N + 41166.9 iV^ - 565.238 - 212.027 A^"' - 16.8356 
-6.53563 7V6 

2.4 200359. + 139453. N + 28750.8 iV^ + 3099.44 - 96.2318 N'^ - 16.3984 

1.5 82989.5 + 55457.5 N + 14428.7 iV^ + 1410.08 + 79.2383 

0,6 13931.7 + 9743.27 N + 2410.02 iV^ + 284.777 A^^ + 12.4061 + 0.284188 
'Yfi 

6.1 -943070. - 161581. iV - 166595. - 9991.42 - 4910.2 iV" + 92.6388 
+ 13.6487 iV^ - 15.9661 iV^ - 5.48745 7V® - 0.997279 TV^ - 0.199456 iV^" 

5.2 -3.12002 10^ - 2.20012 10^ N - 271785. - 102493. + 1777.97 + 334.971 iV^ 
-192.691 - 60.1882 7V^ - 8.77424 - 1.90222 

4.3 -5.8646 10^ - 3.97094 10^ N - 1.15071 10^ iV^ - 22722.2 iV^ - 3636.84 iV* - 814.623 
-258.054 iV^ - 30.149 iV^ - 7.38428 

3.4 -6.00654 10^ - 4.59789 10^ N - 1.11616 10^ iV^ - 144174. + 2260.85 iV* - 257.971 iV^ 
-81.4206 iV^- 17.4821 iV^ 

2.5 -3.73056 10^ - 2.79137 10^ N - 817049. - 94705.5 A^^ „ 6992.55 A^'' + 53.4018 A^^ 
-15.7555 iV^ 

1.6 -1.2573 10^ - 969506. A^ - 282636. A^^ _ 42677.5 A^^ _ 2754.02 A^* - 131.786 A^^ 
0,7 -182134.- 141435. A^- 42673.8 Ar2- 6033.4 Ar3- 464.717 A^*- 14.301 AT^ 

-0.318786 A^6 



Table 4: MZM scheme: coefficients 6^^-, cf. Eq. ( pTTsI ). 
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„w 


4,0,0 


95.0752 









3,1,0 


149.265 


51.9181 


3,0,1 


51.9181 


2,2,0 


128.028 


123.538 


2,1,1 


-9.26743 


2,0,2 


177.083 


-83.8774 






1,3,0 


64.7434 


118.463 


1,2,1 


-20.2218 


1,1,2 


33.8245 


-33.5433 


1,0,3 


33.7500 


0,4,0 


14.6233 


31.5474 


0,3,1 


-2.37189 


0,2,2 


—34.2374 


19.3311 


0,1,3 


OA O ^ P' A 

30.3750 


0,0,4 


— 16.3125 


9.56250 






5,0,0 


-1002.80 









4,1,0 


— 1985.03 


r> A n AAA 

—646.929 


4,0,1 


/-I ^ /-I AAA 

—646.929 


3,2,0 


-2163.00 


— 1771.43 


3,1,1 


-272.688 


3,0,2 


—2308.22 


946.573 






2,3,0 


— 1552.04 


A 1 A A A A 

—2192.02 


2,2,1 


1 A /I '^A P" 

104.725 


2,1,2 


— 1003.60 


A O A A 

608.979 


2,0,3 


P" A A A O A 

—599.238 


1,4,0 


—671.507 


-1 A ^ P" ^ A 

— 1275.70 


1,3,1 


67.695 


1,2,2 


r\r\A m~\c\ 

290.702 


A ^ r" r ci A 

—24.5584 


1,1,3 


—413.455 


1,0,4 


— 166.563 


p' A A rro /I 

59.9734 






0,5,0 


— 121.524 


A <^ <^ A A 

—266.862 


0,4,1 


A t~\ A A A 

8.59444 


0,3,2 


109.555 


-88.119 


0,2,3 


-49.3931 


0,1,4 


-171.978 


46.1564 


0,0,5 


—2.01405 


6,0,0 


13083.1 









5,1,0 


31798.2 


8542.88 


5,0,1 


8542.88 


4,2,0 


42362.6 


27704.8 


4,1,1 


6832.08 


4,0,2 


36182.0 


-13182.8 






3,3,0 


38102.9 


43067.7 


3,2,1 


1007.07 


0,1,2 


26»90.4 


— 10903.4 




12078.4 


2,4,0 


22961.9 


36681.7 


2,3,1 


-686.727 


2,2,2 


1982.93 


-1712.04 


2,1,3 


9048.10 


2,0,4 


9154.83 


-3642.81 






1,5,0 


8272.27 


16001.6 


1,4,1 


-379.749 


1,3,2 


-2399.55 


1194.47 


1,2,3 


1213.86 


1,1,4 


6382.20 


-2878.33 


1,0,5 


91.8928 


0,6,0 


1280.65 


2740.04 


0,5,1 


-52.7711 


0,4,2 


-531.082 


512.566 


0,3,3 


-178.695 


0,2,4 


777.870 


-487.798 


0,1,5 


106.251 


0,0,6 


-156.528 


9.53427 







Table 5: The coefficients c"^;,, c^jj^, and c^;,, cf. Eqs. (| 



,^), O and (I 

corresponding to values of i, j, and k that are not reported are zero. 



The coefficients 
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MS 




MZM 












1,0,0 
0,1,0 






-11/2 
-5/2 






-11/14 
-5/14 


2,0,0 
1,1,0 
0,2,0 


11/16 
5/8 
1/2 


33/8 
7/2 
3 


0.00831444 
0.00755858 
0.00604686 


0.0561224 
0.0510204 
0.0408163 


0,0,2 


3/4 


9/2 


0.00907029 


0.0612245 



3,0,0 
2,1,0 
1,2,0 
1,0,2 
0,3,0 
0,1,2 
"TO" 
3,1,0 
2,2,0 
2,0,2 
1,3,0 
1,1,2 
0,4,0 
0,2,2 
0,0,4 



4,1,0 
3,2,0 
3,0,2 
2,3,0 
2,1,2 
1,4,0 
1,2,2 
1,0,4 
0,5,0 
0,3,2 
0,1,4 

5,1,0 
4,2,0 
4,0,2 
3,3,0 
3,1,2 
2,4,0 
2,2,2 
2,0,4 
1,5,0 
1,3,2 
1,1,4 
0,6,0 
0,4,2 
0,2,4 
0,0,6 



-77/64 
-105/64 
-147/128 

-27/32 
-119/256 

27/64 
0.805664 
1.46484 
7.20703 
13.0078 
8.05664 
1.75781 
2.06909 
-4.13086 
-3.69141 



-6375/128 
-2265/64 
-855/32 
-893/64 
333/32 
273.792 
497.803 
521.554 
381.632 
331.257 
47.2145 
87.3563 
-58.0018 
-14.2917 



0.000692663 

0.000944541 

0.000661179 

0.000485764 

0.00026762 

-0.000242882 

0.00053933 

0.0009806 

0.00119061 

0.00109376 

0.00085394 

0.000150938 

0.000228253 

-0.000209641 

-0.000116359 



-0.0414663 
-0.0565449 
-0.0421311 
-0.0344224 
-0.0152562 
0.00293313 

0.0106715 

0.0194027 

0.027517 

0.0299365 

0.0217858 

0.00437099 

0.00541476 

-0.00713995 

-0.00539322 



-39.5573 
-89.9030 
-129.415 
-99.9105 
-126.293 
-21.1719 
-69.4979 

26.8948 
-9.35655 
-13.9928 

2.66006 
-15.6434 



-2842.37 
-6459.93 
-8195.78 
-4867.47 
-6871.52 
-1702.98 
-3410.39 

622.145 
-728.248 
-698.741 

138.505 
-540.225 



-0.0000621953 

-0.000141353 

-0.000179926 

-0.000107743 

-0.000181477 

0.000105679 

-0.000117862 

0.000126014 

0.00000742886 

-0.0000218022 

-0.00000688957 

-0.0000318276 

0.0000735165 

0.0002005 

0.000404726 

0.000370618 

0.000542556 

0.000278519 

0.000431888 

0.0000308945 

0.0000838964 

0.000182469 

-0.00000319876 

0.000073068 

0.0000308055 

0.0000000587 

0.000020701 

-0.00000211424 



-0.0108866 

-0.0247422 

-0.0319839 

-0.0198859 

-0.0276014 

-0.00595702 

-0.0137125 

0.00503919 

0.000156979 

-0.00263618 

0.000910744 

-0.000466735 

0.00449246 

0.0122522 

0.0224094 

0.0177812 

0.0282151 

0.0115894 

0.0219744 

-0.00073584 

0.00391376 

0.00914719 

-0.000933258 

0.00345242 

0.00149629 

0.00000318554 

0.000740282 

-0.000129075 



Table 6: The coefficients e^-^ and e*^;,, cf. Eq. (O), of the RG functions r]^ and i^t in the 
MS and MZM schemes, to five and six loops respectively. The coefficients corresponding 
to values of i, j, and k that are not reported are zero. 
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2,1,0 
1,2,0 
1,0,2 
0,3,0 
0,1,2 



TF — 



-U.iy08b4 
-0.188964 
-0.100151 
-0.226757 
-0.0306122 
-0.0408163 



TF — 



-0.156841 
-0.219577 
0.122449 
-0.0831444 
0.0589569 



3,1,0 
2,2,0 
2,0,2 
1,3,0 
1,1,2 
0,4,0 
0,2,2 
0,0,4 



0.0837293 
0.135546 
0.123199 
0.162130 
0.0631219 
0.0368576 
0.0128652 
-0.0293008 
-0.0191787 



D 

0.0367151 
0.0941154 

-0.0700545 
0.0972328 

-0.0234689 
0.0271942 
0.0150557 
0.0106478 



2,0 
1,1 

0,2 
0,0 



-0.156841 
-0.0256992 

0.0222978 
0.0181406 



4,1,0 
3,2,0 
3,0,2 
2,3,0 
2,1,2 
1,4,0 
1,2,2 
1,0,4 
0,5,0 
0,3,2 
0,1,4 



-0.043871 
-0.0840225 
-0.0975887 

-0.110406 
-0.0765458 
-0.0421916 
-0.0323891 

0.0254609 
0.00222203 
-0.00529147 
0.00792392 
-0.00144707 



D 

-0.0345052 
-0.0842869 
0.0342725 
-0.100371 
0.0104173 
-0.0600404 
-0.0103217 
0.0047634 
-0.0126929 
-0.00430159 
0.00141906 



3,0 
2,1 

1,2 
1,0 
0,3 
0,1 



0.0367151 
-0.0168042 

-0.0188393 
0.0339000 

-0.000606377 
0.0245772 



5,1,0 
4,2,0 
4,0,2 
3,3,0 
3,1,2 
2,4,0 
2,2,2 
2,0,4 
1,5,0 
1,3,2 
1,1,4 
0,6,0 
0,4,2 
0,2,4 
0,0,6 



0.02V813V — 

0.0698599 

0.101466 

0.0893755 

0.0994876 

0.0600684 

0.0618044 

-0.00398625 
0.0180563 
0.0213326 

-0.00671272 
0.0168769 
0.0030379 

-0.000589662 
0.00331464 

-0.000644392 



D 

0.0131904 
0.0471915 
-0.0273754 
0.0831689 
-0.0164027 
0.0799006 
-0.000641182 
-0.0120957 
0.0375477 
0.000668951 
-0.0124137 
0.00657498 
0.000193903 
-0.00169973 
0.000712328 



4,0 
3,1 

2,2 
2,0 
1,3 
1,1 

0,4 
0,2 
0,0 



-0.0345052 
-0.030022 

-0.0154977 
-0.0337405 
-0.00446734 
-0.0196257 

0.000106648 
0.000711959 
0.00225136 



6,1,0 
5,2,0 
5,0,2 
4,3,0 
4,1,2 
3,4,0 
3,2,2 
3,0,4 
2,5,0 
2,3,2 
2,1,4 
1,6,0 
1,4,2 
1,2,4 
1,0,6 
0,7,0 
0,5,2 
0,3,4 
0,1,6 



-0.0197200 
-0.0560554 
-0.0935780 
-0.0742718 
-0.111549 
-0.068364 
-0.0925957 
-0.0118213 
-0.0354476 
-0.0492935 
0.00447446 
-0.0300979 
-0.0147084 
0.00226375 
-0.00407885 
0.00262247 
-0.00184237 
0.000436103 
0.000962049 
0.00118291 



D 

-0.0147179 
-0.051629 
0.0182722 
-0.0942413 
0.00963234 
-0.107471 
-0.00275039 

0.016807 
-0.0743529 
-0.00164415 
0.0150377 
-0.0274519 
-0.000784818 
0.00239567 
-0.00148025 
-0.00405011 
-0.000306989 
-0.000709295 
-0.00098856 



5,0 
4,1 

3,2 
3,0 
2,3 
2,1 

1,4 
1,2 
1,0 
0,5 
0,3 
0,1 



0.0131904 
0.0038472 

-0.0043311 
0.0318928 

-0.00342079 
0.0166793 

-0.00196111 

-0.00591449 

-0.00290304 

-0.000359204 

-0.00286458 

-0.00126576 



6,0 
5,1 

4,2 
4,0 
3,3 
3,1 

2,4 
2,2 
2,0 
1,5 
1,3 
1,1 

0,6 
0,4 
0,2 
0,0 



-0.0147179 
-0.0226979 

-0.0211958 
-0.0310849 
-0.0132086 
-0.0208751 

-0.00403603 
0.00182088 

-0.00171458 

-0.000383191 
0.00216377 

-0.00290307 

0.0000382403 
-0.0000315805 
-0.000830064 
0.0000197908 



Table 7: The coefficients d'^j^, d'^^jf^, and d^j^, cf. Eqs. (gl^ ), (|4T3|) and 
coefficients corresponding to values of i, j, and k that are not reported are zero. 



The 
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